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Abstract 

Results of a systematic investigation of the vapour-liquid equilibria (VLE) of 30 
individual two-centre Lennard-Jones plus pointquadrupole model fluids (2CLJQ) 
are reported over a range of reduced quadrupolar momentum < Q*^ < 4 and of 
reduced elongation < L* < 0.8. Temperatures investigated are from about 55 % 
to about 95 % of the critical temperature of each fluid. Uniformly the A^jT-l-Test 
Particle Method based on molecular dynamics simulations is used for the generation 
of vapour pressures, saturated densities, and saturated enthalpies. Critical temper- 
atures T* and densities p* are obtained from Guggenheim's equations. Empirical 
correlations for critical data T* and p* as well as for saturated densities p'* , p"* , and 
vapour pressures p* are developed as global functions of the model parameters. In 
most cases they describe the simulation data within their statistical uncertainties. 
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Critical pressures and acentric factors of the 2CLJQ fluid can be calculated from 
present correlations. The present simulations are a sound basis for adjustments of 
the model parameters Q*^, L*, a, and e to experimental VLE data of real fluids. 
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1 Introduction 



Knowledge on vapor-liquid equilibria (VLE) is important in many problems in engineering 
and natural sciences. Among the different ways to model vapor-liquid equilibria, molec- 
ular simulation has the highest potential to yield significant improvements compared to 
existing phenomenological models, especially in terms of predictive power. However, very 
significant further efforts will be needed until molecular simulation based models and tools 
will be sufficiently developed so that their advantages can help process engineers in their 
real world tasks. 

One of the main problems to be overcome is the lack of intermolecular interaction 
models describing vapor-liquid equilibra of real fluids with acceptable accuracy. (With 
"real" fluids we do not refer to non-ideality here, but rather to real world fluids like, 
e.g., 1-propanol or chlorobenzene). As in most applications mixtures are of interest, it is 
necessary to model the different pure fluids with a family of compatible interaction models, 
allowing the application of simple combining rules. Presently the probably most attractive 
and promising family of interaction models are the Lennard- Jones based models, which 
have successfully been applied to a number of real fluids. We mention only some. Vrabec 
et al. have adjusted the parameters of the one centre Lennard- Jones fluid (ICLJ) to 
argon and methane and the two centre Lennard- Jones fluid to ethane [2]. Van Leeuwen 
[3] presented parameters for the one centre Lennard- Jones plus dipole fluid (ICLJD) for 
various real substances. Two centre Lennard- Jones plus dipole (2CLJD) parameters for 
various refrigerants are reported by Kriebel et al. [1] and more recently by Lfsal et al. ^ . 
Finally, the two centre Lennard- Jones plus pointquadrupole (2CL JQ) fluid has successfully 
been used to describe the thermodynamic properties of carbon dioxide (Moller et al. [6J). 
Without reference to real fluids, Dubey et al. [7J have studied phase equilibria of two 
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centre Lennard- Jones dipolar plus quadrupolar fluids (2CLJDQ). O'Shea et al. |8J have 
reported phase equilibria for one centre Lennard- Jones plus quadrupole fluids (ICLJQ). 

The search for an appropriate interaction model for a given fluid is usually a time 
consuming process. In general the focus lies on one pure fluids where the optimization 
of the potential model is done by a number of simulations with subsequent variation of 
the model parameters, cf., e.g., van Leeuwen et al. [9] for methanol. In the present 
work, we follow a new route to develop interaction models that allows fast adjustments 
of model parameters to experimental data for a given class of pure fluids. Here we study 
quadrupolar fluids, of which carbon dioxide undoubted is the most prominent member. 
The favourable results obtained by Moller et al. [6] for carbon dioxide with their 2CLJQ 
model give reason to believe that that model will also be successful for other quadrupolar 
fluids. This is supported by results from a recent study of Calero et al. on refrigerants 
|10] . The idea we follow is to study the thermodynamic properties of the 2CLJQ model 
fluid systematically and in detail over a wide range of model parameters. Using reduced 
coordinates, for the symmetric 2CLJQ fluid only two parameters have to be varied: the di- 
mensionless LJ centre-centre distance L* and the dimensionless pointquadrupole strenght 
Q*^. This allows to cover the entire parameter space of interest with acceptable accuracy 
by studying 30 individual 2CLJQ fluids with different numbers for L* and Q*^. The sim- 
ulation results are correlated empirically in order to be able to interpolate between the 
discrete pairs of L* and Q*^ respectively. 

Based on the results from the present study, it is straightforward to adjust the molec- 
ular interaction parameters of the 2CLJQ fluid to experimental data of real quadrupolar 
fluids. Properties like the critical values of temperature and pressure, the azentric fac- 
tor or the saturated liquid density and vapour pressure are available as functions of the 
molecular interaction parameters. Hence, the development of the molecular interaction 
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model for a given substance is not more difficult than the adjustment of parameters of 
phenomenological thermodynamic models. 

The present work covers the basic molecular simulations for 30 individual 2CL JQ fluids 
and the development of the correlations together with a discussion of those results. The 
application to real fluids will be presented separately. 

Different methods for simulating VLE are available. Many of them are based on Monte 
Carlo techniques (MC) such as the Gibbs Ensemble Monte Carlo by Panagiotopoulos 
[Tl] . that is broadly used. Another method is the Gibbs-Duhem Integration, presented 
by Kofke [121 [13], which uses information from integration of the Clausius-Clapeyron 
equation. Valleau [13] and Kiyohara et al. [15] presented MC based Thermodynamic 
Scaling Methods. Mehta et al. [16] suggested Pseudo-Ensemble Methods that avoid 
particle insertion. Recently, Histogram-Reweighting Methods have been applied for the 
determination of VLE by Kiyohara et al. [TTj [TS] and Conrad et al. [IH] . For simulations 
near the critical points Finite Size Scaling Methods have been applied successfully by 
Wilding [201 EI]. MoUer et al. [22] have suggested the A^pT+Test Particle Method. The 
A^T+Test Particle Method can be based both on MC simulations [23] and on molecular 
dynamics simulations (MD) [6],[2H[25]. The different simulation methods vary in accuracy 
and computing time; we have chosen the A^T+Test Particle Method due to our favourable 
experience with that method especially concerning accuracy. 
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2 Investigated models 



In the present investigation we consider pure two centre Lennard- Jones plus pointquadrupole 
fluids (2CLJQ). It is composed of two identical Lennard- Jones sites a distance L apart 
(2CLJ) plus a pointquadrupole of momentum Q placed in the geometric centre of the 
molecule. The full potential M2Cljq writes as 

wherein 

2 2 
o=l 6=1 

and, as given by Gray et al. [26] 

Uainj, LJi, iV,, ^ 3 Q _ ^ ^ _ ^^^2^2 ^ 2 (s,S,C - 4CiC,)'] , (1) 

with Cfc = cosOk, Sk = sin6k, and c = cos0jj. Herein r^j is the centre-centre distance vector 
of two molecules i and j, Tab is one of the four Lennard- Jones site-site distances; a counts 
the two sites of molecule h counts those of molecule j. The vectors uJi and ojj represent 
the orientations of the two molecules. 9i is the angle between the axis of the molecule 
i and the centre-centre connection line and (pij is the azimutal angle between the axis 
of molecule i and j. The Lennard- Jones parameters a and e represent size and energy 
respectively. 

Among all possible arrangements of the charges in a quadrupole we have chosen the 
four charges to be arranged along the molecular axis in the symmetric sequence +, — , — , 
+ or, having the same energetic effect in pure fluids, —,+,+,—. In the large distance 
approximation the quadrupole interaction reduces to an ideal pointquadrupole interaction 
described by Eq. ([1]). This reduces the number of parameters related to the quadrupole 
to one, namely the quadrupolar momentum Q. 
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Beyond a certain elongation, for small intermolecular distances \rij\ the positive Len- 
nard- Jones term U2clj of the pair potential cannot outweigh the divergence to — oo of 
the quadrupolar term uq. This divergence of M2cljq leads to infinite Boltzmann factors, 
i.e. non-existence of the configurational integral. During molecular dynamics phase space 
sampling within the pressure range in question, this artefact of the 2CL JQ potential causes 
no problem as intermolecular centre-centre distances are very improbable to fall below 
critical values. However, the calculation of the chemical potential by test particle insertion 
often runs into critical intermolecular centre-centre distances. To avoid computational 
problems we followed the proposition of MoUer et al. [6J to put a hard sphere of diameter 
0.4 ■ a directly on the quadrupole site to shield the quadrupolar interaction in critical 
cases. This hard sphere was not active during configuration generation. 

The parameters a and e of the 2CLJQ pair potential were used for the reduction of 
the thermodynamic properties as well as the model parameters L and Q^: T* = Tk/e, 
p* = pa^e, p* = pa^, h* = h/e, L* = L/a, Q*^ = (ea^), At* = At^/^/a. 

The reduced parameters L* and Q*^ were varied in this investigation: L* = 0; 0.2; 0.4; 
0.505; 0.6; 0.8 and Q*"^ = 0; 1; 2; 3; 4. Combining these values leads to a set of 30 model 
fluids to be investigated here. 

In order to achieve a monotonous transition from L* > to L* = we treated 
the spherical fluids ICLJ, where Q*^ = 0, and ICLJQ as two centre LJ fluids with 
L* = 0. This leads to site superposition that is not present when ICLJ and ICLJQ 
fluids are represented by a simple LJ site. Therefore the reduced temperatures, reduced 
pressures, reduced enthalpies and reduced quadrupolar momenta here are fourfold of the 
corresponding values in the one site case. Densities, of course, are not concerned. 

For all simulations the centre-centre cut-off radius Vc was set to 5.0 ■ a. Outside the 



7 



cut-off spliere the ffuid was assumed to have no preferential relative orientations of the 
molecules, i.e., in the calculation of the long range corrections, orientational averaging 
was done with equally weighted relative orientations as proposed by Lustig [27]. The 
quadrupolar interaction needs no long range correction as it disappears by orientational 
averaging. 
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3 Molecular simulation method for VLE data 



For all 30 model fluids the A^T+Test Particle Method (A^T+TP Method) proposed by 
Moller et al. [22l [28] was applied to obtain the VLE data. A detailed description of 
this method has been given by Moller et al. |6j, so we only sketch the key ideas here. 
The ApT+TP Method performs separate NpT simulations in the liquid and the vapour 
phase and uses information about the pressures and the chemical potentials to calculate 
the VLE. Vapour pressures, saturated densities and residual enthalpies h^'^^* {T* , p*) = 
h*{T*,p*) — h^'^* (T*) are evaluated. Test particle insertion for the calculation of the 
residual part of the chemical potential in the NpT ensemble is based on Widom's method 
[29] . Configuration space sampling was done by A^T-molecular dynamics simulations with 
A^ = 864 particles for both liquid and vapour simulations. The dimensionless integration 
time step was set to At* = 0.0015. Starting from a face centred lattice arrangement every 
simulation run was given 10, 000 integration time steps to equilibrate. Data production 
was performed over n = 100, 000 integration time steps. At each production time step 
the information from insertion of 2A^ test particles in the liquid, and A^ test particles in 
the vapour was used to calculate the chemical potential. The dimensionless dynamical 
parameter of ApT-MD-simulations ascribed to the box membrane was set to 2 ■ 10^'^ 
for liquid simulations and to 10^^ for vapour simulations. The high value of A^ allowed 
simulations in the neighborhood of the critical point, and the high value of n was used in 
order to obtain small statistical uncertainties. 

On the basis of the contributions to the free energy A of the 2CLJ interaction, i.e. 
^2CLj from Mecke et al. [30], and of the quadrupolar interaction Aq from Saager et al. 
[3T] . one can construct v42cljq = ^2Clj + ^q- This construction of a hybrid equation 
of state (2CLJQ EOS) is analogous to the method applied to the two centre Lennard- 
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Jones plus pointdipole fluid (2CLJD) by Kriebel et al. [4j. In planning our investigation 
we used this 2CLJQ EOS to estimate the critical temperatures T^*eos (Q*^,-^*) for the 
systems considered. VLE data were determined for temperatures 0.55, 0.60, 0.65, 0.70, 
0.75, 0.80, 0.85, 0.90, 0.925, 0.95 ■ T,*eos(Q*^ ^*)- 

Liquid simulations were performed for the whole temperature range, whereas vapour 
simulations were only performed at temperatures above 0.80 ■ T^f^os- Below this tempera- 
ture the second virial coefficient is sufficient for the VLE calculations. Liquid simulations 
took about eleven hours, vapour simulations about three hours CPU time on a modern 
workstation (e.g. Compaq AlphaStation XPIOOO). 
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4 Simulation results of VLE data and critical data 



Table [T] reports an extract of the VLE data of all 30 model fluids. Given are the vapour 
pressure p*, the saturated liquid density p'*, the saturated vapour density p"*, the residual 
saturated liquid enthalpy h'^^^* , and the residual saturated vapour enthalpy h"'^'^^* for 
temperatures T* ^ 0.55 ■ T*, 0.80 ■ T*, and 0.95 ■ T*. Statistical uncertainties were 
determined with the method of Fincham et al. |32j and the error propagation law. 

Figs. [H [2], [31 H illustrate for Q*^ = 1 and Q*^ = 4 the strong influence of both the 
elongation and the quadrupolar momentum on the 2CLJQ VLE data. Increasing the 
elongation or increasing the quadrupolar momentum strongly influences the shape of the 
density coexistence curve and the slope of the vapour pressure curve. 

When the present 2CLJQ VLE data are compared to the 2CLJQ EOS as described in 
section [HI systematic deviations are observed. This can best be seen in Figs. [2]and[ll For 
models with small quadrupolar momenta and for models with elongations near L* = 0.505 
we find good agreement. This is due to the fact, that the quadrupolar contribution to 
the EOS is based on data of a 2CLJQ model fluid with L* = 0.505. However outside the 
range mentioned above considerable deviations between the EOS and present 2CLJQ VLE 
data are observed. In most cases the EOS underestimates the saturated liquid densities 
and overestimates the saturated vapour densities. Where significant deviations occur, the 
vapour pressures from the 2CLJQ EOS are higher than those from the present study. 

In order to determine the critical data we used the method of Lotfi et al. [24j, who, 
by simple means, found reliable critical data for the ICLJ model fluid. Their values T* = 
1.310 and p* = 0.314 are very close to those indicated later by Potoff et al. [33j . It is known 
that the density-temperature dependence near the critical point is well described by 

1/3 

p* ~ (T* — T*) ' , as given by Guggenheim [31], In order to correlate the saturated 
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densities we followed Lotfi et al. ^24j 

p'* = p* + Ci ■ {T* - T* f^ + C'^ ■ (T; - T*) + C's ■ (T; - T*f/\ (2) 

p"* = p* - Ci ■ (t; - T*)i/3 + C2 ■ (t; - T*) + C3 • (t; - T*f/\ (3) 

The simultaneous fitting of saturated liquid and saturated vapour densities yields not 
only the coefficients Ci, Cg, C3, Cg, C3, but also the critical data p*, T^*. The critical 
temperatures and densities for the 2CLJQ model fiuids are listed in Table [21 In most 
cases present critical temperatures are lower than those estimated by the 2CLJQ EOS. 
Table [2] also contains the critical compressibility factor Zc = p*/ {p*T*). Zc being a non- 
reduced property of the 2CLJQ fluids, it is of particular interest for comparisons to real 
quadrupolar fluids. 

We estimated the uncertainties of T* and p* by refltting VLE data sets that had been 
perturbed stochastically within their own uncertainties. To span the whole range of the 
parameters Q*'^ and L*, we have chosen the model fluids Q*'^ = 0, L* = and Q*'^ = 4, 
L* = 0.8 for these estimations. Ten perturbed VLE data sets for both model fluids have 
been generated. Standard deviations were about cr(Tc) ~ 0.005 and cr(pc) ~ 0.0005. From 
these results we conclude, that the critical temperatures calculated by this method are 
certain up to the second, the critical densities up to the third digit after the decimal point. 
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5 Global correlation of VLE data 

In order to obtain VLE data for the whole range of Q*^, L* and T* we have glob- 
ally correlated the present data. We consider the critical data T*{Q*^,L*), p*{Q*'^,L*), 
the saturated hquid density p'*{Q*'^,L*,T*) and the vapour pressure Inp* (Q*^, L*, T*) 
to be the key VLE data for an adjustment to real fluids. The adequate shape of the 
temperature-density coexistence curve was achieved by simultaneously correlating the 
functions p'* {Q*^,L*,T*) and p"* {Q*^,L*,T*). We want to point out, that it was not in 
the scope of the present investigation to construct a new 2CLJQ EOS. The correlation 
developed here is not designed to compete with an EOS. It shall merely be a working tool 
for a restricted field of application, namely the adjustment of model parameters to data 
of real fluids. The vapour pressure correlation is also used to test the thermodynamic 
consistency of the VLE data from simulations by the means of the Clausius-Clapeyron 
equation. Moreover, the correlations are useful for comparisons with results of other in- 
vestigators. Details of the correlation method are described in the Appendix which also 
contains the resulting correlation functions (cf. Table [3]). 

5.1 Critical properties 

The correlation functions T*{Q*'^, L*) and p*{Q*'^, L*) were assumed to be linear combina- 
tions of elementary functions. The elementary functions and their coefficients are given in 
Table [3] in the Appendix. The quality of the correlations can be studied in Fig. [5l Most 
relative deviations of T* are within 0.5 % which is only slightly more than the individual 
uncertainty estimated before. The critical densities are represented with roughly the same 
quality. It should be mentioned that possible systematic errors may be introduced to the 
critical data by the choice of the exponent 1/3 in the second term of Eqs. ([2]) and ([3]). 
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Both correlations are strictly monotonous in Q*^ and L*. 
5.2 Saturated densities, vapour pressure 

The temperature-density correlation is based on Eqs. and which have five ad- 
justable parameters Ci, C'2, C3, C2, and C3. Another three parameters are introduced 
by the correlation of the vapour pressure (cf. Appendix). Again the functions describing 
the dependency of the correlation parameters on Q*^ and L* were assumed to be linear 
combinations of elementary functions. The elementary functions and their coefficients are 
given in Table in the Appendix. 

A comparison between the correlations and the simulation data can be seen in Figs. [H 
121 El andm A more detailed comparison is given in Fig. El Exemplarily, the relative devi- 
ations of the simulation results from the correlations are shown. Typically, the saturated 
liquid density correlation has the largest deviations for the highest temperature of each 
model, which is due to the large uncertainties in the near critical region. The correlation 
is smooth and shows small deviations in the range of 0.4 % for mid temperatures which 
are most important for adjustments to real fiuids. 

Relative deviations of the saturated vapour densities are not illustrated here. The 
vapour density correlation should not be used below 0.60-Tc(Q*^, L*). At low temperatures 
the correlation is not useful as it does not capture the limiting case of the ideal gas, which 
is independent of the parameters Q*^ and L* . 

In most cases, the vapour pressure correlation represents the simulation data within 
their uncertainties. It has to be mentioned, that the vapour pressures from simulation 
at low temperatures show high uncertainties due to the uncertain values of the chemical 
potential obtained by Widom's test particle insertion in dense liquid phases. 
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By extrapolating the vapour pressure correlation slightly to the critical point, the 
acentric factor [36] 

can be calculated easily from the correlations discussed before, cf. Fig. [71 We find that 
u increases monotonously with increasing Q*"^ and L*. For high Q*^ at approximately 
L* = 0.05 there is a weakly distinct local minimum which probably is an artefact of the 
correlation. Table [2] contains the critical pressures and the acentric factors for the 2CLJQ 
model fluids calculated on the basis of present correlations. 
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6 Discussion 



6.1 Comparison to results of other authors 

Present results are compared with the simulation results of other authors. Fig. [8] presents 
the relative deviations for the saturated liquid density Ap'* = (Pother ~Pcorr)/Pcorr the 
vapour pressure Ap* = (p* other ~ PCTcorr)/PCTcorr- Other authors investigated subsystems 
of the 2CLJQ model fluid, Lotfi et al. the ICLJ [21], Kriebel et al. [2S] and Kronome et 
al. [SZ] the 2CLJ, Stapleton et al. [3H], Smit et al. ^ and O'Shea et al. [8] the ICLJQ, 
and Moller et al. [6] the 2CLJQ. They applied different simulation methods: molecular 
dynamics A^T+TP Method, Monte Carlo A^T+TP Method, or Gibbs Ensemble Monte 
Carlo. 

The saturated liquid densities of the other authors agree in almost all cases within the 
combined uncertainties of our correlation and of the simulation results. The near critical 
points of course show larger deviations. Also for the vapour pressure good agreement is 
observed. 

Present 2CLJQ VLE data show lower statistical uncertainties than those of other 
authors, and they have the advantage that uniformly the same simulation method was 
used to produce them. 

6.2 Thermodynamic consistency test 

The thermodynamic consistency of the simulation data was checked with the Clausius- 
Clapeyron equation 



dlnpl 



AK 



(5) 



dT* p*T*{l/p' 



i/p'*y 
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We used the present vapour pressure correlation in order to evaluate the left hand side 
(LHS) of Eq. ([5]). The right hand side (RHS) of Eq. ([5]) was calculated from our simulation 
data, the uncertainty was calculated by the error propagation law. The requirements of 
Eq. ([5]) are fulfilled at almost all temperatures within the uncertainties of the RHS. We 
conclude that our data are thermodynamically consistent. 

6.3 Locus of the critical point 

The influence of the quadrupole Q*^ and the elongation L* on the critical properties can 
be studied in Table |2] and in Figs. [T] and [21 The critical temperature T* and the critical 
density p* increase together with Q*^ . The increase of T* with Q*^ is more important for 
molecules with small L* than for those with large L* . However, this only holds for the 
absolute increase AT*. The relative increase AT*/T* with increasing Q*"^ only weakly 
depends on L*. This finding is in agreement with results from Garzon et al. [40j, who 
studied the VLE of Kihara fluids with the same pointquadrupole we used. Garzon et 
al. predicted that the 2CLJQ model fluid would qualitatively behave in the same way as 
their Kihara plus pointquadrupole model fluid. 

The critical pressure p* increases with increasing Q*"^ for a fixed elongation L*. The 
critical pressure p* decreases strongly with increasing elongation L* for a fixed quadrupolar 
momentum Q*^. The relative decrease Ap*/p* with increasing L* only weakly depends 
on Q*^. As shown in Figs. [S] and H] the absolute value of the slope of the function Inp* 
vs. 1/T*, i.e. the enthalpy of vaporization Ah*, increases with increasing Q*"^ for a fixed 
elongation L*. It decreases with increasing elongation L* for a fixed Q*^. This is analogous 
to the influence of Q*^ and L* on the locus of the critical pressure p*. It can also be seen 
that the vapour pressure p* of a model fluid of given elongation L* decreases when the 
quadrupolar momentum Q*^ is increased. These results harmonize with those of Garzon 

17 



et al HO]. 



6.4 Deviation from principle of corresponding states 

According to the principle of corresponding states for the 2CLJQ fluid, the curves In (p* /p*) 
vs. T* /T* would form one unique curve regardless of the values of Q*^ and L* . This would 
also be true for the curves T*/T* vs. p* / pi- However, molecular shape and presence of a 
quadrupole do cause deviations from the principle of corresponding states. The widening 
effect of the quadrupole on the density coexistence curve is shown in top of Fig. for 
various 2CLJQ model fluids with L* = 0.6. Molecular anisotropy causes the same effect. 
This is a further confirmation for the conclusion of Garzon et al. [IQ] that VLE coex- 
istence density curves were broadened by the presence of multipoles. The displacement 
of the vapour pressure curves due to the quadrupole is presented in bottom of Fig. 
Qualitatively, the same effect is caused by molecular anisotropy. Fischer et al. [H] have 
reported this effect as well. 

These discernable deviations from the principle of corresponding states are reflected 
by the behaviour of the acentric factor u vs. Q*^ and L*, cf. Eq. (jll). The more a fluid 
deviates from the principle of corresponding states, the higher u will be. 

Garzon et al. [ID], studying the quadrupolar Kihara fluid, introduced the "reduced 
density of quadrupole" 

wherein 

V:,{L*) = (1 + 1.5L*) ■ 7r/6 

is the dimensionless volume of the Kihara molecule, in order to study the increase of the 
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critical temperature 

at; (Qf {Q*\ L*)) = T; (Qf {Q*^ L*)) - T; {Qf (Q*2 = 0, L*)) . (7) 

Within their statistical uncertainties they concluded that AT; vs. forms a unique 
curve for fluids of different anisotropy, i.e. a principle of corresponding states. For the 
2CLJQ fluid we introduce an analogous volume for two intersecting spheres with diameter 
a a distance L* apart 

(L*) = (1 + 0.5 (3T* - L*^)) ■ 77/6. 

Using it in Eq. ([6]) and calculating AT; according to Eq. (JTj), we can not confirm that 
principle of corresponding states for the 2CLJQ fluid, cf. Fig. [TOl top. An alternative 
"effective quadrupolar momentum" suggested by Bohn et al. [12] 

was examined in the same way. The pseudocritical data of unpolar fluids used by Bohn 
et al. [42] have been replaced by present critical data of unpolar fluids. Fig. [TUl bottom, 
shows the apphcation of Eq. ([H]) to Eq. ([7]). Only for high elongations a better accordance 
to the principle of corresponding states can be observed. 
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7 Conclusion 



The present paper is aimed at the quahtative and quantitative improvement of available 
VLE data of the 2CLJQ model fluid. In a systematic investigation the two parameters Q*^ 
and L* of the model fluid have been varied in the ranges to 4 and to 0.8 respectively, 
studying 30 model fluids in detail, including the non-polar and spherical cases. The 
NpT + TP Method has been applied for the production of VLE data in the temperature 
range from about 0.55 ■ T* to 0.95 ■ T* for all 30 model fluids. The comparison of data 
from the present work to a hybrid 2CLJQ EOS has shown the shortcomings of the EOS. 
Single fits of the saturated densities with Guggenheim's equations yielded critical data 
for all 30 systems. 

In order to have useful tools for adjustments of dimensionless VLE data and critical 
data from the present work to real fluids we constructed global correlations of the crit- 
ical data, the saturated density coexistence curve, and the vapour pressure curve. The 
correlations for saturated liquid densities and vapour pressures are sufficiently precise for 
mid-temperatures, which we consider the most important region for the intended adjust- 
ments. Data from the present work show no systematic deviations to results of other 
investigators. In many cases they how ever have lower uncertainties, furthermore the 
whole data set is self-consistent. The Clausius-Clapeyron equation has proven thermody- 
namic consistency of present data. 

Using the correlations from the present work, we studied the influence of Q*^ and L* 
on the locus of the critical point and on the shapes of the saturated density coexistence 
curves and of the vapour pressure curves. A deviation from the principle of corresponding 
states is found, due to the presence of polarity and anisotropy. Our results support those 
of Garzon et al [10]. However, the principle of corresponding states for the increase of 
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the critical temperature of fluids of different anisotropy but equal equivalent quadrupolar 
momentum can not be validated for the 2CLJQ model fluid neither with the "reduced 
density of quadrupole" suggested by Garzon et al. |40| nor with the "effective quadrupolar 
momentum" of Bohn et al. |42j . 

The results from the present work pave the way to applications of molecular simulation 
data to the description of VLE of real quadrupolar fluids. 
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8 List of Symbols 



A coefficient of correlation function 

A free energy 

a interaction site counting index 

B coefficient of correlation function 

b interaction site counting index 

C coefficient of correlation function 

c coefficient of correlation function 

c constant in set of elementary functions 

F function to minimize 

G function to minimize 

h enthalpy 

i data point counting index 

i elementary function counting index 

i particle counting index 

j elementary function counting index 

j particle counting index 

k Boltzmann constant 

k elementary function counting index 

k particle counting index 

L molecular elongation 

LHS left hand side of Clausius-Clapeyron equation 

£ simplified notation for L* 

m mass of particle 

number of particles 

n number of time steps 

p pressure 

Q quadrupolar momentum 

"reduced density of quadrupole" 
"effective quadrupolar momentum" 

q simphfied notation for Q*^ 

RHS right hand side of Clausius-Clapeyron equation 

r site-site distance 

Tc centre-centre cut-off radius 

t time 
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T 


temperature 


u 


inner energy 


u 


pair potential 




molecular volume 


V 


volume 


y 


linear combination or elementary lunctions 


z 


compressibility factor 


a 


coefficient of elementary function 




coefficient of elementary function 


7 


coefficient of elementary function 


A 


relative deviation 


i\CC 


relative deviation of RHb from LHb 


A 7 


enthalpy of vaporisation 




increase of the critical temperature 


At 


integration time step 


d 


statistical uncertainty 


e 


energy parameter 


6 


1 C 11' 

angle of nutation 




elementary function 


P 


density 


a 


size parameter 


a 


standard deviation 




angle of precession 


X 


elementary function 




elementary function 


UJ 


acentric factor 



Vector properties 

r position vector 
u} orientation vector 
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Subscript 



c 


property at critical point 


corr 


from present correlation 




from 2CLJQ hybrid equation of state 


m 


molecular 


other 


from other authors 


Q 


quadrupole 


sim 


from present simulation 


a 


vapour-liquid coexistence 





reference point 


2CLJ 


two centre Lennard- Jones 


2CLJQ 


two centre Lennard- Jones plus quadrupole 



Superscript 

* reduced 

' on bubble line 

" on dew line 

id ideal gas 

res residual property 
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10 Appendix 

10.1 Critical properties 

The functions T*{Q*'^ , L*) and p*{Q*'^,L*) were assumed to be linear combinations of 
elementary functions one of which is a constant c, the others depend either on Q*'^, i.e. 
ipi (Q*^), or on L*, i.e. (L*), or on both, i.e. Xi {Q*^,L*). We restricted the number of 
elementary functions to up to two for both the Q*^- and the L*-dependence and to up to 
four mixed terms. Writing the linear combination of any of the aforementioned functions 
represented by y, we have 

<2 <2 <4 

y (g*^ L*) = c + ^ a, ■ ^,(Q*^) + |^ ■ i^{L*) + |] 7fe ■ Xk{Q*\ L*). (9) 

1=1 j=l k=l 

The crucial step was the selection of elementary functions ipi, ^j, and Xi- The behaviour 
of T* and p* for both the Q*^- and L*-dependence was taken into account for the selection 
of the elementary functions for T*{Q*'^,L*) and p*{Q*'^,L*). The critical data in Table 
12] were fitted to these elementary functions in a usual non-weighted least squares sense. 
The elementary functions and the resulting coefficients are given in Table El 

10.2 Saturated densities, vapour pressure 

The temperature-saturated density correlation is based on Eqs. ((21) and The Q*^, L*- 
dependence was ascribed to T*, p*, and to the coefficients Ci C'^. The correlations for 
T*{Q*'^,L*) and pl{Q*'^,L*) were used in the optimization procedure for the density 
temperature correlations p'*{Q*'^, L*,T*) and p"*{Q*^, L*,T*). The coefficient functions 
Ci(Q*^,L*) ^ C'^i^Q*"^ , L*) were linear combinations of elementary functions in the sense 
of Eq. dni). For the correlation Inp* (Q*^, L*, T*) we restricted the number of coefficient 
functions to three. Several polynomials in 1/T* and an Antoine-type Inp* = A + B/(T* — 
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C) were tested. The best results were obtained for the polynomial ansatz 

1 */^*2 T* rTi*\ /,n*2 T*\ I ^'i(Q ^) ) I '^siQ "^iL) 
lnp^(Q ,L,T)=ci(Q ,L) + h 



which is in accordance to former investigations of Lotfi et al. [2^ and Kriebel et al. 
Again, the coefficients Ci{Q*'^,L*) c^{Q*'^,L*) were assumed to be linear combinations 
of the elementary functions in the sense of Eq. ([9]). 

In principle, in developping the correlations for the saturated densities and the vapour 
pressure, we filtered out elementary functions in the same way as we have done it for the 
correlation of critical data. Fits of VLE data from simulation and from the 2CLJQ EOS 
for some model fluids provided useful information about the topology of the coefficients in 
Q*^- and L*-direction. When an initial set of elementary functions is found, the coefficients 
of these elementary functions can be determined by weighted least squares minimization 
of F and G defined as 

1 



1 

5pp 



mm. 



*^ = ^ x2 [inP'^ L*^T*) - lnp^,sim,i] ^ = min. 

^ (mnp^,i) 

F and G were also used to optimize the initial set of elementary functions. The resulting 
elementary functions and their coefficients are given in Table El 
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List of Tables 



1 Vapour-liquid equilibrium data. Extract from simulation results for 30 
model fluids for low, mid, and high temperatures. At higher temperatures, 
marked by f, data are based completely on simulations. For the remaining 
temperatures the second virial coefficient was used for the vapour phase. 
The numbers in parentheses are the uncertainties of the last decimal digits 
indicated 

2 Critical data, the compressibility factor, and the acentric factor of 30 
2CLJQ model fluids 

3 Elementary functions and their coefficients for the correlations T* (Q*^, L*), 
p*(Q*2,L*), Ci{Q*\L*) to C'^{Q*\L*), and ci{Q*\L*) to c^{Q*\L*). 
The notation is simphfied: £ is L*, g is Q*"^ 
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Table 1: 



rj-i* 




















* 


L* = 0, Q*'^ = 


2.92600 


0.00852 


(33) 


0.82938 


(25) 


0.00299 


(12) 


-26.8247 


(86) 


-0.2197 


(89) 


4.25600t 


0.1469 


(14) 


0.66404 


(55) 


0.04347 


(53) 


-22.457 


(15) 


-2.459 


(33) 


5.05400t 


0.4120 


(34) 


0.4925 


(31) 


0.1475 


(38) 


-17.768 


(77) 


-6.98 


(17) 


L* = 0, g*2 = 1 


2.97000 


0.00796 


(49) 


0.84416 


(22) 


0.00275 


(18) 


-28.4469 


(83) 


-0.212 


(14) 


4.32000t 


0.1472 


(13) 


0.67461 


(58) 


0.04239 


(62) 


-23.448 


(17) 


-2.453 


(43) 


5.13000t 


0.4108 


(34) 


0.4994 


(24) 


0.1398 


(56) 


-18.453 


(60) 


-6.84 


(26) 


L* = 0, g*2 = 2 


3.09100 


0.0057 


(11) 


0.87669 


(24) 


0.00186 


(36) 


-32.4004 


(97) 


-0.165 


(31) 


4.49600t 


0.1417 


(17) 


0.69986 


(61) 


0.03877 


(86) 


-26.051 


(19) 


-2.490 


(77) 


5.33900t 


0.4199 


(36) 


0.5288 


(22) 


0.1375 


(31) 


-20.553 


(59) 


-7.27 


(16) 


L* = 0, g*2 ^ 3 


3.26700 


0.0046 


(11) 


0.91399 


(29) 


0.00144 


(35) 


-37.714 


(14) 


-0.155 


(37) 


4.75200t 


0.1352 


(16) 


0.73335 


(52) 


0.03653 


(67) 


-29.885 


(19) 


-2.745 


(76) 


5.64300t 


0.4213 


(32) 


0.5667 


(17) 


0.1220 


(20) 


-23.694 


(53) 


-7.09 


(11) 


L* = 0, Q*2 ^ 4 


3.47600 


0.0036 


(14) 


0.95228 


(25) 


0.00106 


(43) 


-43.992 


(13) 


-0.144 


(58) 


5.05600t 


0.1202 


(23) 


0.77053 


(46) 


0.0260 


(19) 


-34.733 


(19) 


-2.04 


(27) 


6.00400t 


0.4186 


(42) 


0.6079 


(15) 


0.1082 


(27) 


-27.730 


(54) 


-7.09 


(18) 


L* = 0.2, Q*'^ = 


2.42000 


0.00576 


(31) 


0.73095 


(24) 


0.00244 


(13) 


-22.5710 


(79) 


-0.1737 


(96) 


3.52000t 


0.1092 


(11) 


0.58419 


(52) 


0.03850 


(56) 


-18.792 


(14) 


-2.024 


(47) 


4.18000t 


0.2990 


(22) 


0.4236 


(23) 


0.1279 


(23) 


-14.595 


(52) 


-5.769 


(94) 


L* = 0.2 


, g*^ = 




















2.45850 


0.00514 


(39) 


0.74117 


(23) 


0.00214 


(17) 


-23.7142 


(79) 


-0.158 


(12) 


3.57600t 


0.1092 


(12) 


0.58928 


(50) 


0.03927 


(80) 


-19.412 


(14) 


-2.219 


(69) 


4.24650t 


0.3123 


(23) 


0.4320 


(30) 


0.1317 


(26) 


-15.102 


(69) 


-6.04 


(13) 


L* = 0.2, g*2 = 2 


2.55200 


0.00458 


(49) 


0.76427 


(20) 


0.00183 


(20) 


-26.5697 


(78) 


-0.151 


(17) 


3.71200t 


0.1062 


(12) 


0.60533 


(59) 


0.03580 


(38) 


-21.203 


(17) 


-2.147 


(24) 


4.40800t 


0.3213 


(27) 


0.4448 


(22) 


0.1298 


(27) 


-16.348 


(56) 


-6.21 


(11) 


L* = 0.2 






















2.68950 


0.00300 


(60) 


0.79142 


(26) 


0.00113 


(23) 


-30.474 


(11) 


-0.110 


(22) 


3.91200t 


0.1081 


(16) 


0.62931 


(54) 


0.03498 


(66) 


-23.918 


(18) 


-2.353 


(57) 


4.64550t 


0.3234 


(27) 


0.4559 


(26) 


0.1147 


(41) 


-18.074 


(70) 


-6.11 


(21) 
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Table [TJ continued. 





Pi 




P'* 




p"* 






* 


j^ir res* 




L* = 0.2 


,Q*'' = 


4 


















2.86000 


0.0044 


(15) 


0.81978 


(22) 


0.00158 


(53) 


-35.098 


(11) 


-0.187 


(63) 


4.16000t 


0.1061 


(17) 


0.65261 


(48) 


0.0281 


(14) 


-27.201 


(18) 


-1.94 


(15) 


4.94000t 


0.3428 


(31) 


0.4877 


(20) 


0.1137 


(42) 


-20.837 


(61) 


-6.60 


(24) 


L* = 0.4, Q*'^ = 


1.78750 


0.00325 


(21) 


0.59531 


(20) 


0.00187 


(13) 


-17.2882 


(61) 


-0.1312 


(89) 


2.60000t 


0.06524 


(64) 


0.47244 


(58) 


0.03277 


(80) 


-14.213 


(14) 


-1.708 


(73) 


3.08750t 


0.1892 


(18) 


0.3417 


(26) 


0.1231 


(59) 


-10.955 


(54) 


-5.08 


(21) 


L* = 0.4, Q*2 = I 


1.81500 


0.00267 


(19) 


0.60158 


(20) 


0.00150 


(11) 


-18.0097 


(64) 


-0.1077 


(76) 


2.64000t 


0.06547 


(71) 


0.47589 


(47) 


0.03121 


(59) 


-14.605 


(12) 


-1.626 


(41) 


3.13500t 


0.1910 


(20) 


0.3273 


(46) 


0.1163 


(41) 


-10.845 


(96) 


-4.90 


(16) 


L* = 0.4, g*^ = 2 


1.88100 


0.00226 


(33) 


0.61695 


(19) 


0.00122 


(18) 


-19.8956 


(69) 


-0.096 


(14) 


2.73600t 


0.06747 


(71) 


0.48583 


(52) 


0.03174 


(31) 


-15.731 


(14) 


-1.772 


(22) 


3.24900t 


0.2002 


(20) 


0.3332 


(39) 


0.1166 


(39) 


-11.498 


(91) 


-5.12 


(16) 


L* = 0.4, Q*2 ^ 3 


1.97450 


0.00308 


(45) 


0.63581 


(17) 


0.00160 


(24) 


-22.5113 


(73) 


-0.144 


(22) 


2.87200t 


0.06607 


(81) 


0.49850 


(43) 


0.02904 


(41) 


-17.378 


(13) 


-1.759 


(33) 


3.41050t 


0.2118 


(21) 


0.3515 


(25) 


0.1200 


(36) 


-12.817 


(63) 


-5.69 


(16) 


L* = 0.4, Q*2 ^ 


4 


















2.09000 


0.0023 


(11) 


0.65625 


(20) 


0.00112 


(54) 


-25.654 


(10) 


-0.118 


(57) 


3.04000t 


0.0670 


(11) 


0.51516 


(44) 


0.02688 


(71) 


-19.512 


(15) 


-1.773 


(58) 


3.61000t 


0.2246 


(21) 


0.3635 


(20) 


0.1184 


(34) 


-14.256 


(51) 


-6.07 


(16) 


L* = 0.505, g*2 


= 


















1.55650 


0.00275 


(15) 


0.54062 


(17) 


0.00182 


(10) 


-15.3895 


(52) 


-0.1289 


(72) 


2.26400t 


0.05086 


(49) 


0.42814 


(44) 


0.02815 


(29) 


-12.583 


(11) 


-1.385 


(17) 


2.68850t 


0.1487 


(15) 


0.2927 


(50) 


0.1093 


(38) 


-9.31 


(10) 


-4.42 


(14) 


L* = 0.505, Q*2 


= 1 


















1.57850 


0.00219 


(20) 


0.54667 


(17) 


0.00142 


(14) 


-16.0333 


(55) 


-0.1021 


(97) 


2.29600t 


0.05150 


(62) 


0.43263 


(51) 


0.02783 


(28) 


-12.959 


(13) 


-1.381 


(13) 


2.72650t 


0.1525 


(16) 


0.2972 


(64) 


0.1123 


(40) 


-9.55 


(12) 


-4.59 


(15) 


L* = 0.505, Q*' 


= 2 


















1.63350 


0.00197 


(51) 


0.56070 


(18) 


0.00123 


(33) 


-17.6851 


(69) 


-0.097 


(26) 


2.37600t 


0.05095 


(61) 


0.44022 


(43) 


0.02713 


(70) 


-13.902 


(12) 


-1.470 


(50) 


2.82150t 


0.1611 


(19) 


0.3072 


(31) 


0.1142 


(42) 


-10.273 


(72) 


-4.83 


(15) 
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Table [TJ continued. 





Pi 




o'* 




//* 

y 




^/rcs 


* 


J^ir res* 




L* = 0.505, Q*2 


— 3 


















1.71600 


0.00153 




yj.o I w I 


(^R) 


n nnoQi 

U.UUUc/l 


^99'! 




(70) 
[IV) 


n nsi 

-u.uol 


(90) 


2.49600t 


0.05254 




n 451 84 




0977 


()0) 

[lU) 


-15.321 


OA) 
[11) 


-1 686 




2.96400t 


0.1684 


(21) 


0.3064 


(36) 


0.1138 


(^A) 

[01) 


-10.964 


('901 


-5.17 


(22) 


L* = 0.505, Q*'^ 


— 4 


















1.81500 


0.00089 


V / 




(^R) 

[IS) 




(Of<) 

v / 


99 ^QDS 


[av) 


n n^9 


i''^nl 

OU I 


2.64000t 


0.05430 


(R7\ 
\°' ) 


U.4:U(JU 1 






(77) 


1 7 1 flQ 

± ( . lUU 


Cl f) 

[10) 


1 fi79 


(FlQ) 

Uc* I 


3.13500t 


0.1781 


{lit) 




[01) 


n 1 DQQ 


("^7) 

[01 ) 


-12.274 


(HA) 

[01) 


-5 45 


('181 


L = O.o 


, Q = 




















1.40250 


0.00202 


(^ ^\ 

{10) 


^0090 


n '^) 

[10) 


n nni 4S 


(^ 9) 
[iz) 


±4:. UOOO 


(Af\) 

[10) 


\J. ±U<JZi 


(R9) 

OZ 1 


2.04000t 


0.04288 


(AR\ 
V / 


\J .OcfOOc/ 


[11) 


U .u^ ( 1\J 


(70) 
[ID) 


1 1 4ns 


(^ 0) 

[ID) 


± .oo 1 


UU ( 


2.42250t 


0.1286 


('13'! 


0.2643 


(M) 

[01) 


0.1075 


(33) 


-8.299 


(70) 

[ 1 V) 


-4.26 


fl3l 

[lO) 


L = O.D 


, Q = 




















1.41900 


0.00173 






[I'O) 


n nm 9^ 


(^ 1 


1 4 7n4S 
- 11, i yjio 


('^9) 
[OZ) 




oo ) 


2.06400t 


0.04212 


(A7\ 


u. oy o ±(j 


[It ) 






1 1 7S'i 


(^ 9) 

[IZ) 


1 407 
1 .lyj I 


oo ) 


2.45100t 


0.1301 


(14) 


n 2800 


(20) 


0.1153 


(58) 


-8.804 


(A)) 
[11) 


-4.55 


OH) 

[lO) 


L = O.o 


, Q = 




















1.46850 


0.00136 






07) 


00094 


Of,) 

[ID) 


-1 fi 9369 


('631 


-0 075 

\J .\J i O 


('131 

[10) 


2.13600t 


0.04282 






(A(\) 


U .UOU f 


(^ Q) 
[i\)) 




(^ 'i) 

[10) 




Cl 71 


2.53650t 


0.1311 


(12) 


n 2598 


(65) 


0.1016 


(32) 


-8.79 


OA) 
[11) 


-4.32 


fl2l 


L = O.O 


, Q = 


3 


















1.54000 


0.00144 


(2<^) 


yj .ooo^L 


n 7) 




(^ 7) 


1 S 97'?S 

±0 . ^ 1 oo 


(7A) 
[11) 


u. UOU 


(^ ^) 

[10) 


2.24000t 


0.04270 


(^7) 

[01 ) 


n 41 890 


(A)) 
[11) 


09300 

yj .\j ^fj\J\j 


('33'! 


-1 3 986 


31 

[10) 


-1 989 


('251 


2.66000t 


0.1398 


(l^) 


0.2792 


(AA) 
[11) 


0.1058 


(AR) 

[10) 


-9.87 


fioi 

[IV) 


-4.73 


fl9l 


L* = 0.6 


,Q*'' = 


A 


















1.62800 


0.00083 






^91 
[^1) 




(90) 


9n fi^^ 


(A 1 1 


U. U(JO 


('991 

ZZ 1 


2.36800t 


0.04378 




0.43014 


(AR) 

[10) 


02306 


(70) 

[ 1 U) 


-15.565 


Of) 

[lU) 


-1.495 


(68) 


2.81200t 


0.1487 


(16) 


0.2959 


(20) 


0.1028 


(31) 


-11.108 


(55) 


-4.99 


(14) 


L* = 0.8 


,Q*'' = 




















1.17700 


0.00143 


(15) 


0.43844 


(13) 


0.00124 


(14) 


-12.2840 


(39) 


-0.090 


(10) 


1.71200t 


0.03152 


(36) 


0.34003 


(50) 


0.02282 


(22) 


-9.768 


(12) 


-1.090 


(12) 


2.03300t 


0.0972 


(11) 


0.2103 


(44) 


0.0969 


(32) 


-6.664 


(91) 


-3.71 


(11) 


L* = 0.8 


,Q*'' = 




















1.18800 


0.00120 


(18) 


0.44763 


(13) 


0.00103 


(16) 


-13.0246 


(43) 


-0.078 


(12) 


1.72800t 


0.02961 


(36) 


0.34908 


(40) 


0.02164 


(26) 


-10.278 


(10) 


-1.099 


(12) 


2.05200t 


0.0922 


(10) 


0.2460 


(15) 


0.0874 


(24) 


-7.628 


(38) 


-3.493 


(83) 
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Table [TJ continued. 





Pi 


P 


p'* 


^/rcs 


* 






L* = 0.8 


, Q - 2 














1.22650 


0.00083 (20) 


0.46255 (15) 


0.00069 (17) 


-14.5252 


(60) 


-0.058 


(15) 


1.78400t 


0.02859 (42) 


0.36096 (43) 


0.01992 (37) 


-11.232 


(12) 


-1.109 


(19) 


2.11850t 


0.09230 (95) 


0.2449 (25) 


0.0813 (22) 


-8.067 


(56) 


-3.461 


(87) 


L* = 0.8 


, Q*' = 3 














1.27600 


0.00053 (16) 


0.47977 (14) 


0.00042 (13) 


-16.4793 


(60) 


-0.043 


(13) 


1.85600t 


0.02629 (53) 


0.37721 (36) 


0.0144 (11) 


-12.622 


(12) 


-0.76 


(11) 


2.20400t 


0.0887 (10) 


0.2694 (15) 


0.0699 (34) 


-9.309 


(38) 


-3.29 


(17) 


L* = 0.8 


, Q*' = 4 














1.34200 


0.00033 (26) 


0.49548 (17) 


0.00025 (19) 


-18.6173 


(83) 


-0.030 


(24) 


1.95200t 


0.02440 (67) 


0.39236 (38) 


0.01619 (54) 


-14.226 


(13) 


-1.157 


(71) 


2.31800t 


0.09208 (95) 


0.2893 (14) 


0.0684 (15) 


-10.673 


(36) 


-3.521 


(74) 
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Table 2: 





L* 









0.2 


0.4 


0.505 


0.6 


0.8 









5.236 


4.313 


3.163 


2.735 


2.454 


2.049 


^ C 






0.3143 


0.2740 


0.2251 


0.2032 


0.1850 


0.1577 


Pi 






0.4757 


0.3573 


0.2103 


0.1680 


0.1377 


0.1006 


Pi 






0.2891 


0.3023 


0.2954 


0.3023 


0.3033 


0.3113 








-0.0551 


-0.0194 


0.0319 


0.0439 


0.0553 


0.0731 






1 


5.312 


4.388 


3.195 


2.771 


2.499 


2.097 


C 






0.3169 


0.2796 


0.2259 


0.2036 


0.1863 


0.1619 


Pi 






0.4876 


0.3664 


0.2133 


0.1701 


0.1394 


0.1001 


Pi 






0.2897 


0.2986 


0.2955 


0.3015 


0.2994 


0.2948 


Zc 






-0.0354 


-0.0029 


0.0470 


0.0598 


0.0733 


0.1015 






2 


5.575 


4.554 


3.316 


2.875 


2.580 


2.160 


C 






0.3197 


0.2832 


0.2256 


0.2090 


0.1888 


0.1652 


Pc* 






0.5187 


0.3899 


0.2221 


0.1768 


0.1452 


0.1019 


Pi 






0.2910 


0.3023 


0.2969 


0.2942 


0.2981 


0.2856 


z. 






0.0104 


0.0373 


0.0844 


0.0974 


0.1133 


0.1571 






3 


5.968 


1.807 


3. 189 


3.007 


2.695 


2.266 


C 






0.3303 


0.2880 


0.2309 


0.2097 


0.1940 


0.1680 


Pi 






0.5617 


0.4221 


0.2362 


0.1880 


0.1550 


0.1079 


Pi 






0.2849 


0.3049 


0.2932 


0.2981 


0.2965 


0.2834 


z. 






0.0628 


0.0877 


0.1343 


0.1464 


0.1620 


0.2129 






4 


6.398 


5.143 


3.692 


3.195 


2.866 


2.408 


C 






0.3395 


0.2908 


0.2358 


0.2143 


0.1960 


0.1710 


Pi 






0.5933 


0.4499 


0.2524 


0.2005 


0.1646 


0.1135 


Pi 






0.2731 


0.3008 


0.2899 


0.2928 


0.2930 


0.2756 


z. 






0.1070 


0.1378 


0.1914 


0.2032 


0.2174 


0.2649 
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Table 3: 





Pc* {q. 


c 


1 


0.15075788232-101 


c 


1 


0.31431713778 






0.20472313140 ■ 10"^ 


i^i 


q' 


0.24699993812 - 10^2 




-0.12916709336- 10"^ 




-0.24220105645 • lO'^ 




1/(0.1 + £2) 


0.33194556111 




£2/(0.11 +£2) 


-0.14520352142 


1/(0.1 + £-5) 


0.41364616282 ■ 10"^ 


£7(0-11 + £^) 


-0.42590980830 - lO^^ 


Xi 


q'/iOA+f) 


0.97556493849 ■ 10"^ 


Xi 


fq^iOAl + f) 


-0.27008826493 - 10^2 


gV(0.1 + £5) 


-0.17158397457-10-2 


£5^2/(0.11 + £5) 


0.27854853704 - 10^2 


gV(0.1 + £2) 


-0.81735779939 • IQ-^ 


£2g7(0.11 + £2) 


0.30075661416 • 10"=^ 


gV(0.1+£5) 


0.23012291309 • lO'^ 


£5g7(0.11 + £5) 


-0.50847561388 • 10"=^ 


Ci {q, i) 


C'2 (q, i) 


c 


1 


0.30195487445 


c 


1 


0.49559557035 • 10"^ 


A 




-0.37967456269 • lO'^ 


A 


q' 


0.61257929294 • IO-2 




0.67459200881 • lO'^ 


q' 


-0.16842248337 • IO-2 




£7(^ + 0.4)3 


-0.16082575198 • 10"^ 




e 


0.76651784089 • 10"^ 


£7(^ + 0.4)5 


-0.52279636608 


e 


0.43209147006 • 10"^ 


Xi 


qHy{e + 0.4) 


0.11220680075 • 10"^ 


Xi 


qH'' 


-0.14895101590-10-1 




-0.15237300196 • 10"^ 


<iH' 


-0.21199303602 • 10^2 


qHy{e + OA) 


-0.27116058021 • 10"^ 


qH' 


0.30307914217 • IO-2 




0.39497438901 • lO'^ 








C'^iqJ) 


c 


1 


0.66767184355 • lO'^ 


c 


1 


0.19846809576 • IQ-^ 




q' 


-0.21255206238 • lO'^ 


i^i 


q' 


-0.36770421816 • IO-2 




0.55917439122 • 10"^ 


q3 


0.12514990954 • IO-2 




i 


-0.50078939740 ■ lO'^ 




£2 


-0.59416808329 - lO'^ 


t 


-0.70891584543 ■ 10"^ 


£3 


0.52359791216-10-1 


Xi 


qH 


0.12188553188-10-2 


Xi 


qH' 


0.19662720492- lO^i 


qH^ 


0.86157078089 - 10"^ 


qH^ 


-0.77114540331 - 10^2 


qH" 


-0.14075623914 • lO'^ 


qH' 


-0.26044833098 • IO-2 
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Table [3j continued. 





ci {q,i) 


c 


1 


0.11857649302- 10"^ 


c 


1 


0.43338821199 - 10^ 






0.15564116060- lO'^ 




q' 


0.15036649108 




-0.48721606804- lO'^ 


g3 


-0.20853110369-10-1 






0.42329775623 - 10"^ 




£2/(£2 + 0.75) 


-0.18706065543- 10^ 




-0.28721728138 - lO^^ 


£3/(£3 + 0.75) 


-0.71033866423 


Xi 


qH 


-0.24519442905 - 10"^ 


Xi 


£2g2/(£2 + 0.75) 


-0.57586765382 


qH^ 


-0.52409821736 - 10"^ 


£3g2/(£3 + 0.75) 


0.68025466334 


qH^ 


0.14877168217-10-2 


£2g3/(£2 ^ Q_75) 


0.13582356655 






£3g7(£3 + 0.75) 


-0.19160979995 






c 


1 


-0.26605902479 - 10^ 


c 










-0.11443846237- 10^ 


if^i 


q' 


-0.10592478758 




0.10977797047 




-0.35597308644 - 10-2 




e/{£ + 0.75f 


0.11387292692- 10^ 




£0.5 


-0.88369353198 


£V(£ + 0.75)3 


-0.10829394551 - 10^ 






Xi 


£V/(^ + 0.75)2 


0.11316643469-102 


Xi 






£V/(^ + 0.75)3 


-0.17323579858-102 






£V/(^ + 0.75)2 


-0.16093695721 - 10^ 






£V/(^ + 0.75)3 


0.30266683521 - 10^ 
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List of Figures 



1 Temperature-density coexistence curves of the 2CLJQ model fluids with 
Q*'^ = 1. Symbols: simulation data. Full lines are values from the present 
correlation. Dashed-pointed lines represent the law of rectilinear diameters 
from the present correlation. Dashed lines are values from the 2CLJQ EOS 

[301 [3T] . All errorbars are within symbol size |42 

2 Temperature-density coexistence curves of the 2CLJQ model fluids with 
Q*^ = 4. Symbols: simulation data. Full lines are values from the present 
correlation. Dashed-pointed lines represent the law of rectilinear diameters 
from the present correlation. Dashed lines are values from the 2CLJQ EOS 

[301 EI] ■ AH errorbars are within symbol size |43 

3 Vapour pressure curves of the 2CLJQ model fluids with Q*"^ = 1. Symbols: 
simulation data. Full lines are values from the present correlation. Dashed 
hues are values from the 2CLJQ EOS [SQlEI] |44 

4 Vapour pressure curves of the 2CLJQ model fluids with Q*^ = 4. Symbols: 
simulation data. Full lines are values from the present correlation. Dashed 
hues are values from the 2CLJQ EOS [SOlEI] |45 

5 Relative deviations of simulation data to the correlation. Top: critical 
temperatures AT* = (T*^:^^ — T*^^^^)/T*^^^^. Bottom: critical densities 
^Pc = (Pc.sim - Pc,corr)/Pc,corr- Symbols are: = (•), 1 (■), 2 (♦), 3 

(A),4(v) ke 
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Relative deviations of simulation data to the correlation. Top: saturated 
liquid densities Ap'* = {p'*^^^ — Pcorr)/Pcorr- Full and blank symbols are for 
L* = and L* = 0.505 respectively, with Q*'^ = 1 (•, o), 4 (A, A). Bottom: 
vapour pressures Ap* = {pt,sim~Pa,corr)/Pt,coTT- Full and blank symbols are 
for L* = and L* = 0.8 respectively, with Q*^ = g (•, o), 4 (A, A). . . . 

The acentric factor u of the 2CLJQ fluid vs. Q*^ and L*. Resulting from 
correlation from the present work 

Relative deviations of simulation data of other authors to the correlations 
based on simulations from the present work. Top: saturated liquid densities 

^P'* = (Pother - Pcorr)/Pcorr- Bottom: VapOUr prCSSUrCS p* = (p^^othcr " 

Pa corr)/ Pa corr- Data are taken from: Lotfi et al. p^, ICLJ (o); Kriebel 
et al. [25], 2CLJ with L* = 0.22 (+), 2CLJ with L* = 0.67 (*); Kronome 
et al. [SZ], 2CLJ with L* = 0.67 (x); Stapleton et al. [38j, ICLJQ with 
Q*2 = 4 (4); Smit et al. [33 , ICLJQ with Q*^ = 4 (■); Mdller et al. [B], 
2CLJQ with Q*2 = 3.0255 and L* = 0.65 (•), 2CLJQ with Q*^ = 3.0255 
and L* = 0.779 (A), 2CLJQ with Q*^ = 3.7 and L* = 0.58(T). Errorbars 
show the uncertainties, if they have been indicated 

Deviation from the principle of corresponding states caused by the influence 
of the quadrupole. Top: saturated densities. Bottom: vapour pressures. 
Elongation is L* = 0.6, with Q*^ = (•), 1 (■), 2 (♦), 3 (A), 4 (▼). . . . 
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Deviation from principle of corresponding states for the increase of the 
critical temperature AT^,*, cf. Eq. ((Tj). Top: AT* as a function of the 
"reduced density of quadrupole" and of the elongation L*. Elongations 
areL* = (•), 0.2 (■), 0.4 (♦), 0.6 (A), 0.8 (T). Bottom: AT; as a function 
of the "effective quadrupolar momentum" Q"^ and of the elongation L*. 
Elongations are L* = (•), 0.2 (■), 0.4 (♦), 0.6 (A), 0.8 (T) 



41 



i-i 



Af/% 
c 



Ap /% 
' c 




IF 



A p' 1% 







-1 ■ 



-2 ■ 



-3t 



T 



^ I 



^ ^ A # 



() 



[4 



C) 



1.5 2 2.5 3 3.5 4 4.5 



5.5 6 



20 ■ 



A/7 /% 



10 ■ 
0- 
-10- 
-20- 



() 



o 



o 



4 



^1 — 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



